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The spin Hall effect in a finite ballistic two-dimensional system with Rashba and Dresselhaus spin- 
orbit interaction is studied numerically. We find that the spin Hall conductance is very sensitive to 
the transverse measuring location, the shape and size of the device, and the strength of the spin- 
orbit interaction. Not only the amplitude of spin Hall conductance but also its sign can change. 
This non-universal behavior of the spin Hall effect is essentially different from that of the charge 
Hall effect, in which the Hall voltage is almost invariant with the transverse measuring site and is a 
monotonic function of the strength of the magnetic field. These surprise behavior of the spin Hall 
conductance are attributed to the fact that the eigenstates of the spin Hall system is extended in 
the transverse direction and do not form the edge states. 
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Introduction: The Hall effect (from now on we refer 
as the charge Hall effect) is a well-known important phe- 
nomena in condensed matter physics. It occurs due to 
the Lorentz force that deflects like-charge carriers to- 
wards one edge of the sample creating a voltage trans- 
verse to the direction of current. Recently, another in- 
teresting phenomena, the spin Hall effect (SHE), has 
been discovered and has attracted considerable atten- 
tions. Here, a spin accumulations emerge on the trans- 
verse sides of the sample, when adding a longitudinal 
electric field or bias. If external leads are connected to the 
sides, the pure transverse spin current is generated. The 
SHE can either be extrinsic due to the spin dependent 
scattering^ or intrinsic due to the spin-orbit (SO) inter- 
action. The intrinsic SHE is predicted first by Murakami 
et.al. and Sinova et.al. in a Luttinger SO coupled 3D 
p-doped semiconductor— and a Rashba SO coupled two- 
dimensional electron gas (2DEG)^, respectively. After 
that a number of recent works have focused on this inter- 
esting issue, WWyVVWyVVVQ For an infinite 
system, it was pointed out that the spin Hall conductivity 
is very sensitive to disorder ^Si&Siifliii and the SHE van- 
ishes even in a very weak disorder. On the other hand, 
in the finite mesoscopic ballistic system, the SHE can 
survive>i 2 ii2ii£iiiiii2iiLi& By using the Landauer-Buttiker 
formalism and the tight-binding Hamiltonian22*2i, the 
SHE and spin polarization have been studied in the 
dirtjiiSjiSjM or c i ean i£ii£ii£ mesoscopic samples. These 
investigations show that the SHE is still present below a 
critical disorder. Experimentally, the SHE is observed on 
n-type GaAs 2 ^ and on p-type GaAs2£, where the trans- 
verse spin accumulations are detected by Kerr rotation 
spectroscopy or the circularly polarized light-emitting 
diode, respectively. 

In this paper, we study the nature of the SHE in a finite 
2DEG, and mainly focus on the comparison between the 
SHE and the charge Hall effect. This is because the SHE 
and the charge Hall effect are so analogous, intuitively 



they should have similar properties. In the charge Hall 
effect, the Hall voltage is a universal constant along the 
transverse edge, i.e. it is independent of the transverse 
measuring location and the width of device. At least its 
sign is unchanged. Does the spin current possess similar 
universal behavior in SHE? The results are very surpris- 
ing and show that the transverse spin current in the SHE 
is strongly dependent on the measuring location and the 
device's shaped Not only the intensity but also its sign 
can change. These results indicate that the SHE is not 
as clean as the charge Hall effect, and all the measured 
quantity in SHE are very sensitive on the details of the 
system. We attribute these no-universal behaviors to the 
extensive eigenstates in the transverse direction in the 
SHE. 

model and formulation: The system we considered is 
shown in Fig.l which consists of a finite central ballistic 
region attached to four semi-infinite leads. The Rashba 
and Dresselhaus SO interactions^ are present only in the 
central gray rectangular region with the size N x W. In 
order to study the geometric effect, two zero-SO coupling 
(NSO) zones (central white regions in Fig.l) with the size 
N x m are also patched. All the leads are assumed to 
be clean and ideal, without any SO coupling. The two 
longitudinal leads (lead-1 and lead-2) have the width W, 
which are the same as the width of the central SO re- 
gion. On the other hand, in order to study the local spin 
Hall conductance and its dependence on the measuring 
location, two transverse leads (lead-3 and lead-4) are as- 
sumed to be one-dimensional (ID) with the width 1, and 
they can be coupled to any edge location S along the 
a;— direction. 

The above system can be described by the Hamiltonian 
H = p 2 /2m* + V(x,y) + a(a x p y - a y p x ) + f3(a x p x - 
a y p y ), where a and (3 are the coefficient of the Rashba 
and Dresselhaus SO interactions^ Then in the tight- 
binding representation, this Hamiltonian can be written 
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FIG. 1: Schematic diagram for the mesoscopic four-terminal 
device, which the central gray region (marked by 'SO') has 
the Rashba and Dresselhaus SO interaction, but two central 
white zones (marked by 'NSO') are without the SO coupling. 



FIG. 2: (color online) G s h and G e vs the measuring site S 
for different Vr = 0.03, 0.05, 0.07, and 0.09 with m = [panel 
(a) and (c)], or different width m of the NSO's zone: m — 1, 
3, 5, and 8 with Vr = 0.03 [panel (b) and (d)]. The other 
parameters are Vd = 0, N = 40, and W = 30. 



as 



12,13 



i 

+ E^Mj 

i 

+ H.C 



-t -iV D + Vr 

-iV D - V R -t 

-t -iV R + V D 
-iV R - V D -t 



a i + Syl 



(I) 



where t = h 2 /2m*a 2 is the hopping matrix element with 
the lattice constant a. In order for the band-width of the 
ID lead-3 and lead-4 to be in the same range of —At to 
At, the hopping matrix element in these two leads are set 
to be 2t. Here Vr = ha/2a and Vd = h(3/2a represent 
the strength of the Rashba and Dresselhaus interactions, 
respectively, and Vr and Vd are non-zero only in the 
central gray region. S x and 5 y in Eq.(l) are the unit 
vectors along the x and y directions. 

Since there is no SO interactions in the leads, the spin 
a in the leads are a good quantum number and the def- 
inition of the spin current is unambiguous. Then the 
particle current I pa in the lead-p (p = 1, 2, 3, and 4) 
with spin index a (er =t, or J. stands for the +z or — z 
direction) can be obtained from the Landauer-Buttiker 
formula: I prJ = (e/h) Y. q ^ P ,a> T P<r,qa> (V p - V q )^L where 
V p is the bias in the lead-p and T va ^ qa i is the transmis- 
sion coefficient from the lead-q with spin a 1 to the lead-p 
with spin a. The transmission coefficient can be calcu- 
lated from Tr,, 



pa,qa' — Tr [T pcr G r T qa i G a ] , where the line- 
width function T pa = i(T, prT - 
G r = [G a ]t = {E F -H -j: 
retarded self-energy. After solving JpfT , 



Ep£), the Green's function 
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,„ tJ 1 A and Z r pa is the 
I pai the spin current 
Ip and the charge current I p can be obtained straight- 
forwardly: l e p = e{/ pT + I pl } and Z| = {h/2){I p] - I pl }. 

V and V 2 = 0, 



!p = e{/ P T +/ P i}"and I s p 
The terminal voltages V p are set as: V± 



i.e. a longitudinal bias V is added between the lead- 
1 and the lead-2. The transverse lead-3 and lead-4 act 
as the voltage probes, and their voltages V3 and V4 are 
calculated from the condition /f = Jf = 0. Then the 
transverse spin Hall conductances are: Gs s h — I3/V1 
and GisH — I%jV\. For comparison, we also calcu- 
late the transverse charge currents or the charge con- 
ductances (G^e = I3/V1 and G^e — Zf/Vi) in the same 
device but under a different condition V3 = V4 — in- 
stead of 1^ = 14 = 0. In the numerical calculation, we 
take Ep = — 3.8i which is near the band bottom — At, 
and t = 1 as a energy unit, then the corresponding lat- 
tice constant a w 3nmm& The device's sizes (i.e. N, W, 
and m) are chosen in the same order with the spin pre- 
cession length Lso over the precessing angle ir. Here 
L so = nat/2V R . If taking Vr = OMt, then L so « 50a. 

Numerical results and discussion: First, we consider 
the case that the center region has only Rashba inter- 
action Vr (Vd = 0) and two NSO zones do not exist 
with to = 0. While Vd = 0, it can be shown that 
GzsH — —G^sH = —G s h and Gs e = G± e = G e . The 
spin Hall conductance G s h versus the measuring sites S 
and Vr are depicted in Fig. 2a and Fig. 3a. We see that 
G s h depends on the location of measuring sites S, and 
it can even change its sign, e.g., when Vr = 0.09 (see 
Fig. 2a). On the other hand, at the fixed measuring site 
with different Vr, the curve of G s h versus Vr can also 
cover the range from negative to positive (see Fig. 3a). 
In contrast to the charge Hall effect, their behaviors are 
essentially different. The Hall voltage or the charge Hall 
conductance usually are monotonously increasing func- 
tions of the strength of the external magnetic field. Fur- 
thermore they usually are unchanged with the transverse 
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FIG. 3: (color online) G s h vs Vr for different site S with 
m = [panel (a)], or different width m with 5* = 20 [panel 
(b)]. The other parameters are same with Fig. 2. 



measuring sites. 

Next, we attach two NSO zones to the system (see 
Fig.l). For the charge Hall effect, the charge accumu- 
lates in the transverse boundaries. If two zones having 
no magnetic field are attached, the charge accumulation 
will naturally transfer from the original boundaries to 
the new one, as a result the Hall voltage and the charge 
Hall conductance do not change much. How is the spin 
Hall conductance G s r affected when two NSO zones are 
attached? Fig. 2b and Fig. 3b show, respectively, G s h ver- 
sus the site S and Vr for different thickness m of the NSO 
zone. The results show that the spin Hall conductance 
G s h is strongly affected by the NSO zones. For example, 
in the curves G s h-S, for m = or 1 (see Fig. 2) G s h is 
flat, and it is negative at all site S. With increasing m, 
G s h shows an oscillation behavior. In particular, G s h 
can be positive, i.e. change its sign, for some value of m 
(e.g. m = 5). In the curve of G s h versus Vr it also ex- 
hibits the similar results that G s h is strongly dependent 
on m including changing its sign (see Fig. 3b). 

For comparison, we also show the charge conductance 
G e for the same system but different bias conditions 
V 2 = V 3 = 14 = (see Fig.2c,d). We see that G e is 
always negative and exhibits an oscillation behavior. For 
rn = (i.e. without the NSO zones), G e is weakly depen- 
dent on the site S, whereas for m/0, the oscillatory am- 
plitude of G e increases slightly. In particular, the charge 
conductance G e is nearly ten times larger than the spin 
Hall conductance G s h- This also means G s h is much 
smaller than the universal value Ne/An, where N is the 
channel number, and TV = 2 in the present device because 
that the lead-3's and lead-4's width is 1. Notice that in 
the charge Hall effect, the Hall conductance usually takes 
the universal value Ne 2 /h. 

Let us study the spin Hall conductance G s h versus the 



FIG. 4: (color online) Left panel: G s h (a) and G e (c) vs 
the transverse width W for the site S — 5, fO, and 15. Right 
panel: G s h (b) and G e (d) vs the site S for different transverse 
width W — 44, 45, 46, and 47. The other parameters are 
V R = 0.03, V D = 0, N = 30, and m = 0. 



transverse width W of the center SO's regions. G s h and 
G e versus W exhibit almost periodic peaks (see Fig.4a,c). 
Note that the cutoff energy of the subband (i.e. the 
transverse energy levels) are about n 2 fi, 2 7r 2 /2mW 2 , which 
shifts down with increasing the width W. For the Fermi 
level Ep across a subband, a jump emerges in the curves 
of G s h-W (or G e -W), due to the large density of state 
near the band edge. As a result for a given period (e.g. 
VF=44, 45, 46, and 47), G sH and G e versus the site S 
(see Fig. 4b, d), exhibit the oscillation behavior. As the 
Fermi level across the subband edge (W=47), G s h can 
change its sign while G e are always negative. 

In the following, we investigate the case when the Dres- 
selhaus SO interaction is present, i.e. Vd ^ 0. As men- 
tioned above, at Vd = the spin currents through the 
lead-3 and lead-4 are conserved, i.e., G^sh = —Ga s h- 
However, when Vd ^ and Vr, ^ 0, G^sH ^ —G± s h- 
On the other hand, the spin Hall conductance has the 
symmetry with G^ s h{S) = — G4 S h(W — S) due to the 
symmetry of our system. It is worth to point out when 
Vd = Vr, G4 S h = GisH = 0, which is similar with 
the Ref.(14). In the Fig. 5, Ga s h versus the site S for 
different Vr or different width m of the NSO zone are 
plotted. Here G^h exhibits similar characters as in the 
case of Vd — 0: G^h is very sensitive to the transverse 
measuring site S, and it can even changes its sign (e.g. 
V R = 0.07,0.09). While m ^ 0, G 4sH oscillates with 
the site S along with the variation of its sign. All these 
behaviors are in contrast to the charge Hall conductance. 

Finally, we emphasize that the sensitivity of spin Hall 
conductance to the location of measuring sites is a generic 
feature not due to the ID nature of the lead-3 and lead-4. 
We have performed similar calculations when the width 
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site S (a) 

FIG. 5: (color online) Gash vs the measuring site S for dif- 
ferent Rashba SO coupling strength Vr = 0.03, 0.05, 0.07, 
and 0.09 with m = 0, or different width of the NSO's zone: 
m — 1, 3, 5, and 8 with Vr = 0.03. The other parameters are 
V D = 0.05, N = 40, and W = 30. 

of lead-3 and lead-4 are 3 and 5. The conclusion remains. 
In addition, if the lead-3 and lead-4 are placed at two dif- 
ferent measuring sites along x-direction, G^ s h and Gi s n 



are affected even stronger. 

Why are the characters of the SHE so different with 
the charge Hall effect? Why is the spin Hall conductance 
G s h so sensitive (even its sign) to the measurement site 
5, the shape of device, and so on? We attribute them to 
following two reasons. (1). In the quantum Hall effect 
the edge states emerge and play an important role. How- 
ever for a system that exhibits SHE, e.g., the quasi ID 
quantum wire having Rashba SO interaction, its eigen 
states are extended in the transverse direction and they 
do not form edge states^ (2). The force in the charge 
Hall effect always points to a specific direction, e.g. +y. 
But the force in the SHE is dependent on the spin a, and 
its sign can vary^i 

In summary, the spin Hall conductance is strongly de- 
pendent on the transverse measuring site, the device's 
shape, and the strength of the spin-orbit interaction. Not 
only the magnitude but also its sign can change. These 
characters are very different from that of the charge Hall 
effect, and the spin Hall conductance is not universal as 
the charge Hall conductance. 
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